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Abstract 

The ground state and excitation spectra of the 5 = 1 Heisenberg spin chain with hole 
hopping are investigated by means of the numerical-diagonalization method and by the help 
of the Zhang- Arovas picture. As for the charge sector, two phases are found. It is shown that 
in the weak-hopping region, a phase separation occurs, while in the strong-hopping region, 
the system is in the Tomonaga-Luttinger liquid phase. Irrespective of the presence of the 
spin gap, the present phase diagram is similar to that reported for the 5 = 1/2 t-J chain with 
gapless spin excitation rather than to that for the magnetically-frustrated t-J chain with spin 
gap. For the spin sector, the string correlation is found to develop withstanding the hole 
doping except in a region, whose region is characterized by a generalized string correlation. 
This remains an open question to be studied in future. 
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1 Introduction 



Since the existence of a spin gap on high-T c Cu-oxides was reported [|1], 0, electron systems with 
a spin gap have been studied in the context of the superconductivity. In order to cause a spin 
gap, a number of generalizations of interactions have been made by restricting the dimensionality 
within one dimension: the t-J model with alternating magnetic interaction [0], the t-J model with 
next-nearest-neighbor magnetic frustration [|J, and the double-chain t-J model f|. In particular, 
the double-chain model has lately attracted considerable attention because of the discovery of the 
corresponding materials || 0; see the articles P, H, H3 for a review. 

Amongst such spin-gap states, the Haldane-gap state has been studied extensively. The spin 
gap was conjectured to open for integer-spin antiferromagnetic Heisenberg chains [11, 12|. The 
conjecture was partly confirmed both theoretically |L3| and experimentally [[14]]. Recently, carrier 



doping to the Haldane-gap material was tried ]15j, [U|. The doping causes the reduction of re- 
sistivity and creates states within the Haldane gap ]17|- There are several theoretical proposals 
to describe such experimental results 



19, 2G 



These proposals are based upon some model 
Hamiltonians in which each hole is regarded to have the inner freedom of the spin S = 1/2, and 
propagates along an S = 1 spin chain. The hopping amplitude of carriers is, however, estimated 
to be so small that only the considerations upon the impurity problems could explicate essential 



properties 21 



In the present paper, we investigate the following one-dimensional S — 1 t-J model under the 



periodic-boundary condition: 



W = -fE E (iW + h-C-) + ^ESrS I+ i. 

i=l o-=-l,0,l 



i=l 



The operators and &f CT obey the hard-core-boson statics 

0. 

Sera', 



The operators {Sj} are expressed as 




(2) 



(3) 



where the matrices s = (s x , s y ,s z ) represent the 5 = 1 spin operators. This model may be too 
simplified to describe the above magnetic material. Apart from the practical interest, however, 
the present Hamiltonian ([[]) may be useful for studying essential effects of a spin gap to phase 
diagrams. This is quite accordant with our motivation explained at the beginning of the present 
paper. 

For some limiting cases, some properties of the ground state of the Hamiltonian (P are known 
as follows. At the filling n — 1, the model reduces to the 5 = 1 spin system; the system has the 
so-called Haldane gap, AE 



0.4105 J |26[. The two-particle problem corresponding to the 
low-density limit, namely n — j,, can be solved exactly; two particles make a bound state in the 
region J/t > 1, while they are deconfined in the region J/t < 1. Hence, it is expected that in the 
region J/t > 1 particles are apt to form an island of particles for arbitrary filling n; the so-called 
phase separation takes place. Hence, the characteristics of the system is rather simple in this 
phase-separated region. Notice that this situation is relevant to the above material. We stress 
here that the present paper mainly concerns the liquid phase J/t<l, in which the properties of 
the fluid are nontrivial. 

The present paper is organized as follows. In the next section, we explain the Zhang- Arovas 
argument |E7| which is also applicable to the present problem ([I]). This yields a rough feature 
of the phase diagram in our system, and predicts the existence of a spin gap withstanding the 
hole doping, and so forth. Referring to these predictions, we present in §3 numerical results for 
the model ([!]) both on the charge and spin sectors. It is shown that such a new phase as is not 
predicted from the Zhang- Arovas argument spreads in a low-density region. In addition to this, 
we report that the present system in the metallic phase is not essentially affected by the presence 
of the Haldane gap. The phase diagram resembles that of the 5=1/2 t-J model |28] without 
spin gap rather than that of the magnetically-frustrated t-J -J' model || with a spin gap. The 
present spin-liquid state seems to be rather unique compared to other spin-liquid states. In the 
last section, we give a summary of the present paper. 



2 Zhang- Arovas picture for the present model 



In this section, we explain the Zhang and Arovas picture relevant to the present system | 27| . They 
discussed the following system; 



H 



-*E E (Cw + h.cO + jEfs.-s.+x + ^s.-s^) 2 

i=l <r=-l,0,l i=\ V 6 



(4) 



The magnetic interaction including the biquadratic term is precisely the same as that in the AKLT 



model [g^, Q , which can be solved exactly at the ground state 

First, we consider the impurity problem t = 0. At the condition t = 0, the model reduces to 
a lot of AKLT sectors under the open-boundary condition; the system is cut off at the positions 
of holes to form independent finite AKLT spin chains. On the other hand, the ground state of 
the AKLT model is called the valence-bond solid (VBS) state. The VBS state for the system of 
length I under the open-boundary condition is given as [[H]] 

|VBS> = (alH&t) 1 -* f n(aj6} +1 - &]a] +1 ) j (ajy(bj)^\0), (0 < p, q < 1), (5) 
where the S = 1 spin operators are expressed in terms of the Schwinger representation 

= ajbi, S~ = Oib\, SI = i 1 - — — and afaj + b\bi = 2. (6) 

The formula ([5]) can be interpreted such that a valence bond is formed over each bond in a 
systematic manner. The edge freedom of p and q in eq. (|5|) can be expressed in terms of the 
following different representation. The edges holding hole(s) between them can be classified in 
terms of the total magnetization of these edges. The magnetization S 2 = implies that a singlet 
valence bond is formed over the hole(s), while the magnetization S 2 = 1(1 + 1) implies that a 
triplet bond is formed; see Fig. 1. Once the hopping amplitude t is assumed to be nonvanishing, 
this expression becomes essential. Hence, it is consequently found that there are two types of 
holes that can be classified with respect to the inherent bond; we call the former (latter) 'singlet 
hole' ('triplet hole'), hereafter. 

Second, we consider the system (|j) with a finite hopping term t. It was reported that the 
singlet-hole state with the wave number k is an exact eigenstate, while the triplet-hole state 



is not [p7 |. Hence, the Hamiltonian for the singlet holes can be given exactly in terms of the 



spinless-fermion operators; 

2 

n = -tj2 (4 c i+i + h - c -) - o J 4ci4 +1 c i+ i (7) 

i i 

Hence, the band of one singlet hole is given by = — 2t cos k. That for one triplet hole is 
approximately estimated as = y cosfc f27|. Hence, the singlet-hole state is concluded to be the 



ground state in the presence of the hopping term. The Hamiltonian (|7|) can be transformed in 
terms of the Jordan- Wigner transformation into the following ferromagnetic Heisenberg chain: 



J 



W = -4f £ ( a*a? +1 + afaf +l + -a* a 



'i+i 



where the Hilbert space is restricted within J2i °f = 2N—L; N denotes the number of the particles. 
Now, we are in a position to grasp a rough feature of the phase diagram of the present system 
with singlet holes. In the XY limit J/t — > 0, the ground state is in the XY phase, i.e., the 
Tomonaga-Luttinger liquid phase. In the Ising limit J/t — > oo, the state is phase-separated: the 
positive-magnetization sector and the negative-magnetization sector extend separately. Along the 
conditions J2i °f — =F 2, i.e., the high- or low-density limits, the state shows the confinement- 
deconfinement transition at J/t = 3. Note that this transition point coincides the isotropic 
condition of the spin system (|8|). Keeping in mind of the particle-hole symmetry q <-> (— l)*cj of 
the Hamiltonian (0), the phase diagram of the singlet-hole state may be drawn schematically as 
in Fig. 2. 

Lastly, we explain the properties of the spin sector. If the singlet-hole state is realized at the 
ground state, the string correlation |E 



O z stlins U ~i) = (S*^=* S *S*) (9) 



develops withstanding the hole doping. The reason is as follows. We assume that the singlet-hole 
state should be an ground state. Then, focusing only on the S = 1 particles, i.e., disregarding the 
holes, the spin sector is expressed in terms of the VBS state; a valence bond is formed over each 
neighboring 5=1 particles. On the other hand, the string correlation (|9]) is known to develop on 
the VBS state. 

The existence of the developing string correlation can be interpreted as the existence of a spin 
gap because of the following reason. This existence implies that the state should be relevant to 
the VBS state. A finite energy is required to excite the state, whose energy corresponds to exiting 
one of the valence bonds to a triplet bond. Hence, in the present paper, we study the string 
correlation rather than the spin gap directly. 



3 Numerical Results and Discussions 

In this section, we investigate the ground-state properties of the system (ffl) by means of the 
exact-diagonalization method. The results are discussed by the help of the Zhang- Arovas picture 
explained in the previous section. Phase diagrams are given both for the charge and spin sectors. 



3.1 Elementary excitations 

In this subsection, we show dispersion relations both for the phase-separated region J/t>l and 
liquid region J/t<l, before estimating the precise phase boundary. Dynamical structure factors 
are evaluated so that they can yield the respective elementary excitations of the charge and 
spin sectors. This separation enables us to analyze the charge sector in terms of the Tomonaga- 
Luttinger theory; see the article pl|| for a review. 



First, in Fig. 3 (a), we present the dispersion relation for the system with L = 14, J/t = 10 ^> 1 
and n = 13/14 which is a single-hole system in the phase-separated region. The low-lying bands 
originate in the single-hole dispersions. Because of the degeneracy of the band with respect to 
the total magnetization see Fig- 3 (a), the wider (narrower) band is seen to correspond 

to singlet (triplet) one. According to the Zhang-Arovas picture explained in the previous section, 
the band width of the singlet hole is approximately three-times wider than that of the triplet 
one. The dispersions in Fig. 3 (a) actually agree with this prediction. As for the spin sector, we 
find that the spin gap opens independently of the wave number k, AE ~ J. The reason of this 
independency is due to the condition J/t <C 1, where a hole can hardly propagate and can be 
regarded as an impurity. As a consequence, the state is independent of the wave number. 

Second, in Fig. 3 (b), we show the dispersion relation for the system with L = 14, J/t = 1 and 
n = 12/14 which is a system with a hopping amplitude comparable to the magnetic interaction. 
The structure of elementary excitations is no more manifest. The spin and charge elementary 
excitations are separated with use of the following analysis. 

We employed the dynamical structure factor S(k, u) both for the spin and charge sectors, 

.Wgc(A^) = J2( \ n -k\j)(j\n k \0)S{co - Ej + E ), (10) 

3 

and 

S spiR (k,uj) = J2(0Sl k \j)(j\S z k \0}5(uj - Ej + E ), (11) 

3 

where n k = J2i,a e ~ lkTi b\ a h,a and S% = J2i,o e ~ lkTi Si- The factors show 5-function peaks 
at the levels which have the above particular transition probabilities. They are evaluated with 



use of Mori's continued-fraction-expansion formalism |35], pq| . The factors (|10D and ( JlT|) for the 
system with L = 14, J/t = 1 and n = 12/14 are shown in Fig. 4 (a) and (b), respectively, where 
the peaks are broadened with the parameter r\ = 0.04 such as S(x) = — ilm^-. As was discussed 



in the previous section, the charge sector in the region J/t<l is expected to be described in terms 
of the Tomonaga-Luttinger theory. It is crucial to estimate the sound velocity v c that is defined 
with respect to the low-lying dispersion relation uj = v c k [k ~ 0). Through collating Fig. 4 (a) 
with Fig. 3 (a), it can be seen that the sound velocity for the charge elementary excitation should 
be estimated with respect to the level A; v = [Ex — E )/ (x)- The estimated sound velocity for 
the system with L — 14, n — 12/14 and J/t varied is shown in Fig. 5. 

As for the spin elementary excitation depicted in Fig. 4 (b), the first exited state locates at 
the wave number k = k-p = mr. This location coincides with that for the S = | t-J model, whose 
model was, however, shown to have a gapless spin excitation spectrum by the use of the rigorous 
solution along the tractable line J/t = 2 ||37|| . The spin excitation of the present model is shown 



to have a spin gap in the subsection 3.3. 

3.2 Phase diagram of the charge sector 

In this subsection, we investigate the phase diagram of the charge sector. It is analysed in term 
of the Tomonaga-Luttinger liquid theory, where the estimation of the sound velocity v c given in 
the previous subsection plays an essential role. 

The exponent K p which governs the Tomonaga-Luttinger liquid state can be expressed by 
§§, 

Kp = l Vc n 2 K, (12) 

where k denotes the compressibility k = ^^§§r- This compressibility is estimated by means of 
the difference method as 

L E [N + 2) + E [N-2)-2E [N) 
K =W 4 • (13) 

Every correlation exponents of the charge sector are governed in terms of the single parameter K p ; 
for example, the charge- density wave and the pairing correlations obey the following power-law 

(OLsit y (OOdensit y (0)) oc l/r K ?, (14) 

<OU(r)O pair (0)> oc 1/r 1 ^, (15) 

(OdensityO') = Ea Aa h ia Opair(«) = Ea4a4+l,a 

respectively. At K p = 1, the above correlations ( |14| ) and (pl|) decay with the same exponent. The 
former (latter) should be dominant, if the liquid is "repulsive (attractive)." Hence, it is seen that 
the magnitude of the exponent K p indicates a degree of the attraction among the particles and 
that the particles are "noninteractive" at K p = 1. This interpretation will be found to be valid 
through referring to the exact estimation of K p for various models |4l], |4"T|. 



In Fig. 6, we presented the estimation of l/K p for the system with L — 14, n — 12/14 and 
J/t varied. In the weak-interaction limit J/t — > 0, only concerning the charge sector, the present 
system converges to the spinless-hard-core gas. For such a system, the exponent K p is given by 
K p = 0.5 exactly. This is quite consistent with the present numerical estimation in Fig. 6. As 
the magnetic interaction J/t is increased, the exponent K p increases. In the region l<J/t<1.6, 

the exponent is estimated as K p > 1. Hence, the liquid is attractive; the paring correlation (|T5|) is 
dominant in this region. On the other hand, the exponent K p is negative in the region J/t> 1.6. 

Therefore, the compressibility k is also negative in this region; see eq. (0)- Hence, the gas 
is unstable; we can regard the system as being phase-separated in this region. Note that these 
features are consistent with the argument explained in §2. 

With use of the similar data for various values of the filling n, we obtained the phase diagram 
for the charge sector, see Fig. 7. As for the charge sector, two qualitatively- different phase 



diagrams have been reported so far (see Fig. 8 (a) and (b)). The phase diagram of the present 
model belongs to the type depicted in Fig. 8 (a). Notice that the phase diagram depicted in 
Fig. 2 which is given according to the Zhang- Arovas argument also belongs to the type of Fig. 8 
(a). The phase diagrams (a) and (b) were reported to hold for the one-dimensional S — 1/2 t-J 
model p8f and for the t-J model with the next-nearest-neighbor magnetic interaction J' = 0.5J 



[|J, respectively. The former model shows a gapless spin excitation. On the other hand, the latter 
model is expected to have a spin gap near the filling n ~ 1, because the model converges to the 
Majamder-Ghosh model in the limit n — > 1. In the parameter region with n ~ 1 and J/t < 1, 



the diagram Fig. 8 (a) shows that the liquid has the exponent K p = 0.5, while the diagram (b) 
shows that the liquid is either phase-separated or has the exponent K p ^> 1. Hence, the liquid of 
the latter case is strongly attractive one. The reason has attributed to the presence of a spin gap, 
so far. 

The present model, however, belongs to the type (a) in spite of the presence of a spin gap, 
called the Haldane gap. It suggests that the RVB pattern is crucial for enhancing the pairing 
correlation; the present RVB pattern belongs to the type called the VBS state (see the previous 
section). 

3.3 Spin sector 

In this subsection, we concentrate on the spin sector of the present model ([]]). The development 
of the string correlation (§) indicates the presence of a short-range RVB state which is essentially 
the same as the VBS state; the state is thus concluded to have a spin gap corresponding to exiting 
one of the valence bonds to a triplet. Physical properties of the liquid phase are of greater interest, 
because the phase-separated particle system is no more than an island of the S = 1 spin chain in 
a sea of holes. 

In Fig. 9, we have plotted the string correlation 0* tring (L/2) for the systems with J/t = 1.0 
and L, n varied. As is expected, the correlation decreases with the filling n decreased. The 
plot, however, does not show clearly whether it develops or not. In order to see it definitely, we 
calculated the Binder parameter |43|, H5| for the string corelation; 



U(L) = 1-3^, (16) 

O = i^^Sti^. 
L i 

With the system size increased, the Binder parameter increases if the corresponding correlation 
is developing. We show the plot for the systems with J/t = 0.2, 1.2, 2.2 and L, n varied (see Fig. 
10 (a), (b) and (c), respectively). At J/t = 0.2, the string correlation develops withstanding the 
doping. Holes coherently move in the VBS state without disturbing it. Hence, at J/t = 0.2 we 
can conclude that the singlet-hole picture is valid at the ground state. The above behaviour of the 
spin gap is in contrast with that for the frustrated S = 1/2 t-J chain M. For this model, it was 

3 

reported that the spin gap closes in an over-doped region n<— . At J/t = 1.2, in the low-density 

region n ~ 0, the correlation is found to be of short range. This phase is absent in the Zhang- 
Arovas picture introduced in the previous section. The nature of this new phase is discussed later. 
The correlation is seen to develop at J/t = 2.2, which is rather strong. It is quite natural, if we 
notice that the system is phase-separated (see the phase diagram depicted in Fig. 7). 

With use of similar data for various values of J/t, we draw schematically the transition line 
in Fig. 11. As is discussed in §1, the two-particle system can be solved exactly; there exists a 
confinement-deconfinement transition at J/t = 1. We have utilized this additional fact for drawing 
the transition line in Fig. 11. 



In order to see the nature of the above new phase, we generalize the correlation ([|) to the 
following form [46], 

0^(\i-M = (5?e"£E (17) 
This has been successfully applied for the spin-5 antiferromagnetic Heisenberg chains [46|, ^7|, [IB] 
and the 5 = 1 spin chain with the alternating interaction |3E|, . 

For the 5 = 1 spin chain with the alternating interaction, the following was reported. The 
correlation at 9 = n is maximal in the Haldane phase, while the correlation at 9 = tt is vanishing 
and that around 9 = tt/2 is maximal in the dimer phase. In Figs. 12 (a) and (b), the correlation 
O z (L/2,9) is plotted for the systems with (a) L = 14, J/t = 0.5 and n = 12/14 and (b) L = 14, 
J/t = 1.2 and n = 4/14, respectively. In Fig. 12 (a), a maximum locates at 9 = tt, while in Fig. 
12 (b), a maximum locates about at 9 = tt/2. The former actually confirms that the state is in 
the Haldane phase. On the other hand, the latter implies that the state is, somehow, in the dimer 
phase; the 5 = 1 particles may make a singlet bound states in pair. A similar state is reported 
in the low-density region for the 5 = 1/2 t-J chain [OSJ, |50| . In the present model, however, the 
region spreads over a large region. It is noteworthy that the state has a spin gap corresponding 
to exiting one of the singlet bound pairs to a triplet pair. As a consequence, all the phases that 
appear in Fig. 12 are found to show a spin gap. 



4 Summary 

The ground state of the one-dimensional 5 = 1 t-J model ([I]) is investigated. In the weak- 
magnetic-interaction region J/t < 1.6, the ground state is in the liquid phase, while the system 

is phase-separated in the strong-interaction region J/t > 1.6. In spite of the presence of the spin 
gap, the phase diagram of the charge sector belongs to the type that was reported for the one- 
dimensional 5=1/2 t-J model rather than that for the frustrated t-J model with a spin gap. 
Although the presence of a spin gap has been speculated as a sign of strong attraction among the 
particles, the present result suggests that the pattern of the RVB formation is rather crucial. As 
for the spin sector in the liquid region, in addition to the phase conjectured by Zhang and Arovas, 
a new phase is found to appear. The phase spreads in the low-density region, and the string order 
is disturbed. The disturbed state, however, possesses the generalized string correlation of 9 = |. 
It suggests that the 5=1 particles make bound states in pair. A similar phase is observed in the 

5 = 1/2 t-J model, in whose model, however, it expands in a very limiting region [T2^, |5D| . It turns 
out consequently that a spin gap outlives to open in the whole region. The present spin-liquid 
state seems to be rather unique compared to other spin-liquid state. 
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Figure captions 

Fig. 1: Schematic drawing of the generalized VBS state in the presence of holes. Over the 
hole(s), either a singlet or a triplet bond is supposed to be formed. 

Fig. 2: Schematic drawing of the phase diagram of the Hamiltonian (f|) due to the Zhang- 
Arovas argument. 

Fig. 3 (a) (b): Dispersion relations for the systems with (a) L — 14, J/t — 10 and n = 13/14, 
and (b) L — 14, J/t — 1 and n = 12/14, respectively. The symbols +, x and o denote the levels 
of the quantum number J2i S£ = 0, ±1 and ±2, respectively. 

Fig. 4 (a) (b): Dynamical structure factor S(k,u) for (a) charge excitation and (b) spin 
excitation for the system L = 14, J/t = 1 and 12/14 in Fig. 3 (b). 

Fig. 5: Velocity of the charge-excitation spectrum for the system L = 14, n = 12/14 and J/t 
varied. 

Fig. 6: Inverse of the exponent K p of the charge sector for the system L = 14, n = 12/14 and 
J/t varied. 

Fig. 7: Phase diagram of the charge sector. The contours of the exponent K p are also shown. 
The plots +, x and o denote the data estimated for the systems of the length L = 10, 12 and 14, 
respectively. 

Fig. 8: Qualitatively different phase diagrams of the charge sector reported so far. The type 
(a) was reported for the S = 1/2 t-J model with gap less spin excitation E?|, while the type (b) 
was reported for the frustrated t-J- J' model with spin gap 

Fig. 9: String correlation Og tring (L/2) for the system with J/t = 1, and n varied. 

Fig. 10 (a) (b) (c): Binder parameter U(L) of the string correlation for the systems with (a) 
J/t = 0.2, (b) J/t = 1.2 and (c) J/t = 2.2, and n varied. 

Fig. 11: Schematic drawing of the phase diagram of the spin sector. 

Fig. 12 (a) (b): Generalized string correlation Cg tring (L/2, 6) against the angle 6. It is calcu- 
lated for the systems with (a) L = 14, J/t = 0.5 and n = 12/14, and (b) L = 14, J/t = 1.2 and 
n = 4/14. 
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